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Summary 


The  objective  of  this  study  is  determination  of  the  accuracy 
with  which  the  station-source  backazimuth  can  be  determined  from 
single  station,  three-component  data  in  the  case  of  industrial 
explosions  at  near  regional  distance.  The  single  station  backazimuth, 
derived  from  the  polarization  of  the  P-wave,  was  a  parameter 
reported  by  U.S.  stations  during  the  GSETT-2  experiment.  One  of 
these  stations  (BLA),  located  at  Blacksburg,  Virginia,  is  situated  to 
record  numerous  large  mining  explosions  in  the  distance  range  100 
to  300  km.  The  explosions  can  be  reliably  located  using  arrival  time 
data  from  a  regional  network. 

A  moving  window  polarization  analysis  of  the  explosion  P- 
wave  arrivals  was  performed.  In  a  fashion  similar  to  the  procedures 
followed  during  the  GSETT-2  experiment,  maximum  rectilinearity  of 
motion  was  used  as  a  criterion  for  judging  the  "best”  backazimuth 
estimate  from  the  three-component  P-wave  data.  The  resulting  best 
estimate  backazimuth  is  compared  with  that  derived  from  the 
independently  located  mines. 

The  results  of  the  comparisons  show  that  three-component 
signal/noise  ratios  less  than  2.0  can  result  in  very  unreliable 
backazimuths.  However,  for  S/N  ratios  exceeding  2,0,  the  mean 
backazimuth  error  is  6  degrees  with  standard  deviation  21  degrees. 
The  error  in  the  determinations  does  not  appear  to  be  strongly 
correlated  with  signal/noise  ratios  in  the  range  2  to  10. 
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Characteristics  of  the  signals  which  contribute  to  the  errors  in 
the  single-station  backazimuth  estimates  include; 

1.  Very  emergent  initial  motions  from  the  delay-fired 
explosions. 

2.  Off-azimuth  P-wave  arrivals  and/or  converted  phases 
arriving  very  early  in  the  P-wave  coda. 

3.  Steep  apparent  angles  of  incidence,  averaging  22  degrees 
from  vertical,  which  tend  to  reduce  the  signal/noise  ratios 
on  the  horizontal  components. 

Optimum  P-wave  data  segments  for  reliable  source 
backazimuth  estimates  are  restricted  to  very  short  time  intervals  (1 
second  or  less)  beginning  with  the  initial  P-wave  motion. 

A  subset  of  data  featuring  Rg  arrivals  was  analyzed  using  a 
phase  difference  criterion  to  estimate  the  station-source 
backazimuth.  The  mean  backazimuth  error  for  27  events  is  -2 
degrees:  the  standard  deviation  is  26  degrees,  which  is  only  S 

degrees  larger  than  the  result  obtained  from  the  P-wave  arrivals. 
This  suggests  that  Rg  can  be  used  to  reliably  estimate  the  source 
direction  in  some  situations. 
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1. 


Introduction 


As  part  of  the  GSETT-2  experiment,  a  high-quality  six -channel 
(three  component  short-period  and  broadband)  seismic  system  was 
installed  at  Blacksburg,  Virginia,  in  January,  1990.  The  Blacksburg 
location  is  well  situated  to  record  industrial  explosions  in  the 
distance  range  100  to  300  km  arising  from  surface  coal  mining 
operations  in  Ohio,  West  Virginia  and  Kentucky. 

This  report  examines  the  accuracy  with  which  the  source- 
station  backazimuth  can  be  determined  at  station  BLA  from 
polarization  analysis  of  the  mining  explosion  generated  P-wave  and 
Rg-wave  signals.  As  a  basis  for  the  accuracy  assessment,  the 
backazimuths  derived  from  the  three  component  analysis  at  BLA  are 
compared  with  the  locations  of  the  mines  estimated  using  arrival 
time  data  recorded  at  stations  of  the  Virginia  Tech  and  Tennessee 
Valley  Authority  regional  seismic  networks. 

2.  The  Blacksburg  Three<Component  Station 

The  Blacksburg,  Virginia,  GSETT-2  station  was  installed  in 
January,  1990  in  the  WWSSN  vault  on  the  campus  of  Virginia  Tech. 
This  vault  is  founded  on  Cambrian  dolomite  of  the  Appalachian 
Valley  and  Ridge  geologic  province.  The  station  is  equipped  with 
Teledyne  Geotech  GS-13  short-period  and  BB-13  broadband  sensors. 
Digitization,  multiplexing  and  calibration  are  performed  by  a 
Teledyne-Geotech  RDAS-200  unit,  and  time  synchronization  is 
obtained  via  a  Kinemetrics,  Inc.  Omega  receiver/clock.  The  short- 
period  and  broadband  sample  rates  are  40  and  10  samples/sec, 
respectively.  The  data  are  demultiplexed,  event-detected  and 
archived  by  a  Science  Horizons,  Inc.  NOMAD  workstation. 

2.1  Calibration 

Relative  differences  in  system  response  on  the  two  horizontal 
components  result  in  systematic  errors  in  backazimuth  estimates. 
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Hence,  the  system  responses  of  the  short-period  channels  were 
determined. 

The  RDAS-200  unit  provides  for  the  capability  of  either  pulse- 
type  transient  or  sinusoidal  steady  state  calibration  test  signals  to  be 
applied  to  the  sensor  calibration  circuits.  The  transient  calibration 
pulses  were  used  to  determine  the  short-period  system  response 
using  the  approach  described  by  Chapman  et  al.  (1988).  Essentially, 
the  method  fits  an  ideal  damped  oscillator  amplitude  response  to  the 
Fourier  amplitude  spectrum  of  a  recorded  system  calibration  pulse. 
A  nonlinear  least  squares  algorithm  is  used  to  estimate  the  system 
period,  damping  and  gain. 

The  amplitude  responses  of  the  horizontal  channels  were  found 
to  differ  by  less  than  2  percent  in  the  1  to  10  Hz  frequency  band. 
The  vertical  component  response  differs  from  that  of  the  two 
horizontals  by  approximately  9  percent.  These  differences  were 
removed  by  mathematical  correction  of  the  recorded  data. 

3.  The  Regional  Network 

The  network  stations  used  to  determine  the  locations  of  mining 
explosions  are  shown  in  Figure  1.  The  six  stations  in  western 
Virginia  and  in  West  Virginia  are  operated  by  Virginia  Tech,  whereas 
the  stations  in  Tennessee  and  Kentucky  are  part  of  a  network 

operated  by  the  Tennessee  Valley  Authority  (TV A).  The  Virginia 

Tech  stations  are  recorded  digitally  in  an  event  triggered  mode.  The 

TV  A  stations  are  recorded  on  16  mm  Develocorder  photographic  film. 

4.  Location  and  Characteristics  of  the  Mining  Explosions 

Figure  1  shows  the  locations  of  the  37  mining  explosions  used 

in  the  study.  Examination  of  the  signals  using  trace  overlays 

recorded  by  the  Virginia  Tech  network  indicates  that  the  data  set 

consists  of  multiple  explosions  at  nine  mine  sites  plus  single 
explosions  at  nine  additional  sites.  Confirmation  of  the  geographic 


locations  of  the  explosions  through  contact  with  the  mine  operators 
has  been  achieved  for  three  mines.  These  three  confirmed  locations 
are  the  Central  Ohio  Mine  in  southeastern  Ohio,  the  Ruffner  Mine  in 
West  Virginia  and  the  Starfire  Mine  in  eastern  Kentucky,  shown  in 
Figure  2. 

The  epicenters  of  the  explosions  were  located  using  the 
combined  Virginia  Tech  and  TV  A  network  arrival  time  data  (both  P 
and  S‘Lg  arrival  times).  The  location  program  HYPOELLIPSE  (Lahr, 
1980)  was  used  with  a  three  layer  crustal  velocity  model  (Bollinger 
et  al.,  1980).  The  network  locations  and  location  quality  estimates 
are  given  in  Table  1.  In  cases  where  the  event  epicenter  was  not 
known  from  independent  information  and  where  multiple  events 
were  recorded  from  the  same  location,  the  network  location  with 
smallest  statistical  uncertainty  was  adopted  as  the  epicenter  of  the 
group  of  explosions. 


Figure  1.  Network  stations  (including  three-component  station  BLA) 
are  shown  by  circles.  Crosses  indicate  mine  locations 
yielding  multiple  explosion  signals.  Mines  for  which  only 
a  single  explosion  was  recorded  are  shown  by  triangles. 


3 


Table  1 


EVENT 

NO. 

LAT 
(Deg.  N1 

LONG  ERH2 
(Deg,  m  (kml 

DIST 

fkm^ 

BAZ 

(Deg.l 

BAZ 

(EsO 

S/N 

A 

39.80 

81.69 

308 

339 

334 

0.88 

1.1 

21 

A 

39.80 

81.69 

- 

308 

339 

368 

0.91 

1.4 

31 

A 

39.80 

81.69 

- 

308 

339 

20 

0.93 

1.7 

35 

A 

39.80 

81.69 

- 

308 

339 

239 

0.87 

1.7 

36 

A 

39.80 

81.69 

- 

308 

339 

331 

0.89 

2.2 

10 

B 

37.43 

83.12 

- 

241 

277 

275 

0.99 

3.0 

13 

B 

37.43 

83.12 

- 

241 

277 

250 

0.98 

7.2 

16 

B 

37.43 

83.12 

- 

241 

277 

270 

0.96 

7.2 

28 

B 

37.43 

83.12 

- 

241 

277 

289 

0.97 

4.1 

44 

B 

37.43 

83.12 

- 

241 

277 

243 

0.97 

3.8 

15 

C 

81.81 

- 

139 

298 

248 

0.94 

6.1 

47 

C 

81.81 

- 

139 

298 

271 

0.98 

3.8 

19 

- 

38.53 

80.39 

15 

146 

1±4.2 

357 

0.88 

1.8 

22 

- 

38.53 

80.39 

1 1 

146 

1±4.2 

9 

0.85 

3.2 

12 

- 

38.13 

82.40 

7 

202 

301±1.5 

299 

0.98 

6.9 

39 

- 

38.13 

82.40 

7 

202 

301±1.5 

286 

0.98 

5.0 

43 

- 

38.13 

82.40 

21 

202 

301±1.5 

313 

0.94 

3.1 

1  8  -  37.99  82.54  5  206  295±0.7  280  0.97  3.3 

3  2  -  37.99  82.54  3  206  295±0.7  265  0.98  9.5 


24 

- 

37.58 

82.18 

8 

161 

285±2.8 

144 

0.93 

1.4 

38 

- 

37.58 

82.18 

9 

161 

285±2.8 

283 

0.93 

1.3 

41 

- 

37.58 

82.18 

10 

161 

285±2.8 

314 

0.96 

5.4 

Table  1  (continued) 


EVENT 

NO. 

Minel 

LAT 
(Deg.  N) 

LONG 
(Deg.  W) 

ERH2 

(km) 

DIST 

(km) 

BAZ 

(Deg.) 

BAZ 

(Est) 

17 

38.03 

81.69 

9 

144 

309±1.6 

307 

0.96 

4.4 

42 

- 

38.03 

81.69 

9 

144 

309±1.6 

311 

0.91 

2.5 

34 

- 

38.03 

81.69 

5 

144 

309±1.6 

320 

0.90 

2.5 

27 

- 

81.88 

14 

154 

303±2.1 

287 

0.94 

3.0 

33 

- 

81.88 

7 

154 

303±2.1 

308 

0.97 

8.7 

8 

- 

38.34 

82.78 

17 

243 

302±3.6 

308 

0.92 

2.0 

9 

- 

37.52 

82.64 

99 

200 

280±2.3 

279 

0.96 

4.8 

1  1 

- 

38.04 

82.26 

1 1 

187 

300±2.2 

273 

0.95 

3.8 

14 

- 

37.94 

82.65 

6 

212 

293±1.1 

264 

0.96 

6.9 

25 

- 

37.19 

81.51 

12 

97 

269±4.5 

274 

0.94 

3.7 

26 

- 

37.85 

82.54 

7 

200 

292±1.7 

284 

0.90 

3.3 

29 

- 

38.28 

82.36 

4 

208 

305±1.0 

279 

0.97 

7.9 

37 

- 

38.28 

81.21 

6 

137 

330±2.2 

336 

0.89 

2.4 

40 

- 

38.04 

82.16 

19 

179 

304±3.7 

252 

0.97 

3.2 

=  Central  Ohio  Coal  Company; 

B  =  "Starfire"  Mine,  eastern  Kentucky; 

C=  "Ruffner"  Mine,  Arch  Minerals,  Inc.,  Yolyn,  West  Virginia; 

-  =  Mine  location  derived  from  network  arrival  time  data;  not  confirmed 
by  independent  information 

^ERH  =  Horizontal  error  measure  from  location  program  HYPOELLIPSE. 

^Errors  estimated  on  the  basis  of  the  angles  subtended  by  the  HYPOELLIPSE 
68%  confidence  error  ellipses. 
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The  location  capabilities  of  the  network  were  examined  by 
comparing  the  network  derived  epicenter  locations  of  the  Starfire, 
Ruffner  and  Central  Ohio  Mines  with  their  known  locations.  Figure  2 
shows  the  HYPOELLIPSE  68%  chi-square  confidence  regions  for  the 
horizontal  projections  of  the  hypocenter  locations.  Also  shown  in  the 
figure  are  the  actual  locations  of  the  three  mines.  The  network 
locations  for  the  explosions  at  the  Starfire  Mine  in  Kentucky  and  the 
explosions  at  the  Ruffner  Mine  in  West  Virginia  are  in  good 
agreement  with  the  actual  locations.  The  network  location  capability 
for  the  Ohio  explosions  is  degraded  because  of  very  low  signal-to- 
noise  ratios  for  that  more  distant  source. 

Figure  3  shows  the  location  error  ellipses  for  the  remaining 
events  in  the  data  set.  With  the  exception  of  Event  9  in  Kentucky, 
which  was  well  recorded  only  by  the  Virginia  Tech  network,  the 
error  ellipse  semi-major  axes  are  ail  less  than  25  km.  The 
uncertainty  in  the  network  derived  backazimuths  must  be  taken  into 
consideration  in  order  to  assess  the  quality  of  the  single-station 
estimates.  The  horizontal  angles  subtended  by  the  HYPOELLIPSE 
error  ellipses  as  viewed  from  station  BLA  were  calculated  and  are 
listed  in  Table  1.  The  maximum  uncertainty,  as  defined  by  the 
angular  extent  of  the  error  ellipse,  is  ±4.6  degrees,  for  Event  #25. 
This  error  estimate  does  not  account  for  any  systematic  bias  due  to 
inaccuracy  of  the  velocity  model.  However,  this  latter  source  of  error 
appears  to  be  small,  for  as  shown  in  Figure  2,  the  calculated  locations 
of  the  Starfire  and  Ruffner  Mines  are  near  the  true  mine  locations. 

Information  obtained  from  the  mine  operators  at  the  three 
mines  indicates  that  the  larger  explosions  usually  are  the  result  of 
"cast  blasting,”  a  technique  which  typically  utilizes  delay-fired, 
multi-row  charge  designs,  similar  to  the  explosions  studied  by 
Chapman  et  al.  (1992). 
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Figure  2.  Confirmed  mine  locations  are  shown  by  the  circles. 

Location  of  station  BLA  shown  by  the  triangle.  Ellipses 
indicate  the  horizontal  projections  of  the  68%  chi-square 
confidence  ellipsoids  for  blast  epicenters  from  location 
program  HYPOELLIPSE. 
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Figure  3.  Confirmed  mine  locations  are  shown  by  th^  circles  with 
associated  letters  as  listed  in  Table  1.  The  location  of  the 
three-component  station  is  shown  by  the  triangle.  The 
ellipses  indicate  the  horizontal  projections  of  the 
HYPOELLIPSE  68%  chi-square  confidence  ellipsoids  for 
mine  locations  inferred  from  blast  epicenters. 


5.  Method  of  P*Wave  Polarization  Analysis 

The  P-wave  backazimuth  estimates  from  the  three-component 
BLA  station  were  derived  using  the  method  developed  by  Jurkevics 
(1988).  The  processing  is  carried  out  in  the  time  domain  using  a 
series  of  overlapping  time  windows.  For  each  time  window,  a  3  x  3 
data  covariance  matrix  is  formed  from  the  auto-  and  cross-variances 
of  the  three  components  of  motion.  The  principle  axes  of  the 
polarization  ellipsoid,  (A.i  Ui,  i  =  1,  2,  3)  for  each  window  are 
determined  by  solving  the  eigenproblem  for  eigenvalues  (Xi,X2and 
X3)  and  eigenvectors  (ui,  U2,  and  U3).  For  P-wave  motion,  the 
azimuth  of  propagation  is  estimated  from  the  horizontal  orientation 
of  the  eigenvector  corresponding  to  the  largest  eigenvalue,  ui: 

Pazimuth  =  tan'*  (U2i/u3i),  (1) 

where  uji,  j  =1,2,3  are  the  direction  cosines  of  the  eigenvector 
associated  with  the  largest  eigenvalue.  A  measure  of  the  degree  of 
rectilinearity  of  the  particle  motion  is  given  by 

R=  1 (2) 

which  is  1.0  when  X2  and  X3  are  zero,  for  pure  P-wave  motion.  The 
method  can  be  made  frequency  dependent  by  applying  the  analysis 
to  narrow  bandpassed  time  series  data. 
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6.  P’Wave  Data  Analysis 


The  BLA  short-period  data  were  corrected  for  minor 
differences  in  system  response  of  the  three  component  channels. 
Frequency  domain  assessment  of  the  P-wave  to  pre-P-wave 
signal/noise  levels  was  made  by  calculating  the  Fourier  amplitude 
spectrum  of  the  P-wave  and  pre-P-wave  data  segments  using  10  sec 
time  windows.  Explosions  with  signal/noise  ratios  less  than  2  in  the 
frequency  range  1  to  10  Hz  were  rejected.  The  three  component  data 
were  then  bandpass  filtered  using  a  four-pole  Butterworth  filter  with 
comer  frequencies  1.0  and  10.0  Hz. 

The  duration  of  the  moving  time  window  used  in  the 
polarization  analysis  is  a  critical  parameter  for  determining  accurate 
backazimuth  estimates  from  the  P-wave  arrivals.  Its  value  will 
depend  upon  the  bandwidth  of  the  signal,  the  sample  rate  and  the 
polarization  characteristics  of  the  signal.  Trial  and  error  testing 
showed  that  a  window  duration  of  O.S  seconds  generally  provided 
the  best  backazimuth  determination  with  our  data  set.  Given  the 
data  bandwidth  of  1  to  10  Hz  and  a  sample  rate  of  40 
samples/second,  windows  less  that  0.5  seconds  in  duration  gave 
unstable  results  due  to  the  small  number  of  sample  points,  whereas 
larger  windows  tended  to  reduce  resolution  because  of  the  non¬ 
stationary  polarization  character  of  the  coda  immediately  following 
the  initial  P-wave  arrival. 

The  polarization  analysis  was  performed  in  an  automated 
manner,  with  the  best-estimate  of  the  backazimuth  determined  on 
the  basis  of  maximum  rectilinearity  near  the  time  of  initial  P-wave 
motion. 
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Examples  of  the  various  steps  involved  in  the  processing 
procedure  are  shown  in  Figures  4  through  6.  Shown  in  Figure  4  are 
the  data  from  Explosion  #16  at  the  Starfire  Mine,  rotated  into 
vertical,  radial  and  transverse  components.  Figure  5  shows  the 
Fourier  amplitude  spectra  of  the  10  second  P-wave  and  pre-P-wave 
noise  segments.  Figure  6  shows  a  5  second  duration  sequence  of  the 
data  centered  near  the  time  of  P-wave  onset,  with  continuous 
running  estimates  of  backazimuth,  rectilinearity  and  signal/noise 
ratio.  The  signal/noise  quantity  was  calculated  as  the  ratio  of  the 
average  three-component  amplitude  in  the  S  second  moving  window 
to  the  average  long-term  three-component  amplitude  using  a  10 
second  data  segment  prior  to  P-wave  onset. 


I  *■  '  '  'T  '  r  '  ^  r*  i‘'  r  '  •  '  i*  '  i"  f 
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Figure  4.  Time  series  data  for  Explosion  #16  at  the  Starfire  Mine  in 
Kentucky.  The  data  have  been  rotated  into  vertical  (top), 
radial  (center),  and  transverse  (bottom)  components. 
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Figure  5.  Fourier  amplitude  spectra  of  the  vertical  component  P  wave 
amval.  Explosion  #16  (solid  line),  and  the  pre-P  wave  noise 
(dashed).  Ten  second  time  windows  were  used. 
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Figure  6.  Results  of  polarization  processing  for  Explosion  #16.  From  the 
top,  the  traces  are:  (1)  Vertical,  (2)  radial,  (3)  transverse,  (4) 
station-source  backazimuth  (deg),  (5)  rectilinearity,  and  (6) 
three-component  signal/noise  ratio.  Note  that  for  rotated 
components,  the  true  backazimuth  angle  is  180  degrees. 

Vertical  lines  indicate  time  at  which  the  rectilinearity  is 
maximum.  The  duration  of  the  moving  window  for  polarization 
analysis  is  O.S  seconds. 


13 


7.  Results  of  the  P>Wave  Analysis 

The  backazimuth  errors,  defined  as  the  differences  in  the 
network  derived  estimates  and  the  results  from  the  three- 
component  P-wave  polarization  measurements,  are  plotted  as 
functions  of  signal/noise  ratio  and  rectilinearity  in  Figures  7  and  8, 
res];)ectively.  Figure  7  suggests  that  ratios  less  than  2  can,  but  do  not 
always,  result  in  very  unstable  estimates  of  the  polarization 
parameters,  as  evidenced  by  the  two  points  with  backazimuth  errors 
in  excess  of  90  degrees.  These  two  events  are  #35  and  #24  (Table  1). 
Note  that  other  examples  of  explosions  at  both  of  these  mines, 
featuring  similar  signal/noise  ratios,  produced  backazimuth 
estimates  much  closer  to  the  true  values.  An  example  involves 
Events  35  and  20,  both  at  the  Central  Ohio  Mine.  Figures  9  and  10 
show  the  three-component  recordings  for  the  explosions.  Figures  11 
and  12  show  the  vertical  component  P-wave  and  pre-P-wave  noise 
spectra  of  the  two  signals,  and  Figures  13  and  14  summarize  the 
three -component  processing  results  for  both  events.  Note  that  the 
moving  window  backazimuth  estimates  near  the  times  of  P-wave 
onset  are  stable  at  values  near  the  true  value  of  180  degrees  for  the 
rotated  components,  in  the  case  of  Event  20  (Figure  14),  whereas 
those  for  Event  35  (Figure  13)  change  rapidly  at  times  near  the  P- 
wave  onset.  Clearly,  the  random  changes  in  the  character  of  the 
noise  may  severely  affect  the  backazimuth  determination,  when  the 
signal/noise  ratio  is  small. 

Large  signal/noise  ratios  do  not  guarantee  accurate 
backazimuth  estimates,  at  least  for  BLA  and  this  mine  explosion  data 
set.  Figure  7  shows  no  evidence  for  the  backazimuth  error  to 
decrease  with  increasing  signal/noise  ratio,  in  the  range  2  to  10. 
Deleting  the  two  observations  with  errors  in  excess  of  100  degrees 
results  in  a  mean  azimuth  error  of  -6  degrees;  the  standard  deviation 
is  21  degrees. 


14 


Backozimuth  Error  (Deg)  ^  Bockozimuth  Error  (Deg) 


Three-component  station  backazimuth  error  versus 
signal/noise  ratio. 


Rectilinearity 


Figure  8.  Three-component  station  backazimuth  error  versus 
rectilinearity. 
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Figure  9.  Time  series  data  for  Explosions  #35  at  the  Central  Ohio 
Mine.  (Top)  Vertical  component,  (middle)  radial 
component,  (bottom)  transverse  component. 


l»fe«  UK 


ac  sa 

Figure  10.  Hme  series  data  for  Explosion  #20  at  the  Central  Ohio 
Mine.  (Top)  Vertical  component,  (middle)  radial 
component,  (bottom)  transverse  component. 
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Figure  11.  Fourier  amplitude  spectrum.  Explosion  #35:  Vertical 
component  P  wave  arrival  (solid)  and  pre-P  wave  noise 
(dashed).  Ten  second  time  windows  were  used. 
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Figure  12.  Fourier  amplitude  spectrum.  Explosion  #20:  Vertical 
component  P  wave  arrival  (solid)  and  pre-P  wave  noise 
(dashed).  Ten  second  time  windows  were  used. 
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Figure  13.  Results  of  polarization  processing  for  Explosion  #35.  From  the 
top,  the  traces  are:  (1)  Vertical,  (2)  radial,  (3)  transverse,  (4) 
station-source  backazimuth  (deg),  (5)  rectilinearity,  and  (6) 
three-component  signal/noise  ratio.  Note  that  for  rotated 
components,  the  true  backazimuth  angle  is  180  degrees. 

Vertical  lines  indicate  time  at  which  the  rectilinearity  is 
maximum.  The  duration  of  the  moving  window  for  polarization 
analysis  is  0.5  seconds. 
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Figure  14.  Results  of  polarization  processing  for  Explosion  #20.  From  the 
top,  the  traces  are:  (1)  Vertical,  (2)  radial,  (3)  transverse,  (4) 
station-source  backazimuth  (deg),  (S)  rectilinearity,  and  (6) 
three-component  signal/noise  ratio.  Note  that  for  rotated 
components,  the  true  backazimuth  angle  is  180  degrees. 

Vertical  lines  indicate  time  at  which  the  rectilinearity  is 
maximum.  The  duration  of  the  moving  window  for  polarization 
analysis  is  0.5  seconds. 
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Figure  8  shows  the  backazimuth  error  plotted  as  a  function  of 

rectilinearity.  Again,  there  is  no  tendency  in  this  data  set  for  the 

errors  to  decrease  with  increasing  values  of  this  parameter  in  the 
range  0.8  to  1.0.  Thus,  it  appears  that  the  P-wave  arrivals  from 
some  of  the  mine  explosions  are  arriving  from  off-azimuth  directions, 
possibly  due  to  refraction  from  laterally  heterogeneous  velocity 
structure. 

As  a  rule,  the  explosion  data  set  exhibits  very  emergent  and 
complex  P-wave  arrivals.  The  emergent  character  of  the  signal  may 
be  due  to  the  extended  nature  of  the  delay-fired  source-time 
function.  For  example,  typical  large  explosions  at  the  Starfire  Mine 
involve  up  to  60  separate  charges,  and  the  duration  of  the  source 
time  function  can  be  in  excess  of  1  second  (Chapman,  et  al.,  1992). 

This  is  also  the  case  for  some  of  the  larger  explosions  at  the  Ruffner 

Mine  (personal  communication,  mine  engineers)  and  is  probably 
typical  of  most  of  the  explosions  in  our  data  set  from  other  mines  as 
well.  Results  from  the  polarization  analysis  indicate  that,  in  general, 
only  the  initial  portion  of  the  P-wave  arrival  exhibits  strongly 
rectilinear  motion,  which  is  polarized  in  the  source-station  azimuth. 
Off-azimuth  arrivals  consisting  of  converted  and  scattered  energy 
appear  at  a  very  early  stage  in  the  coda:  rarely  can  accurate 
estimates  of  the  P-wave  backazimuth  be  obtained  from  portions  of 
the  coda  at  times  greater  than  1.0  second  following  the  initial  onset 
time.  Also,  as  shown  in  Figure  15,  the  P-wave  angles  of  incidence 
with  the  vertical  are  small,  averaging  22  degrees  with  a  standard 
deviation  of  9  degrees.  This  results  in  a  reduced  signal/noise  ratio 
on  the  horizontal  components,  which  is  approximately  one-half  of  the 
total  three-component  value  given  in  Table  1.  This  phenomenon, 
which  may  be  station  dependent,  tends  to  further  reduce  the 
accuracy  of  the  backazimuth  determination  for  low  S/N  events. 
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Incidence  Angle  (Deg) 


Figure  15.  Apparent  angle  of  incident  P  wave  versus  three- 
component  signal/noise  ratio. 

Figures  16  through  27  show  five  examples  with  signal/noise 
ratios  exceeding  4  wherein  the  polarization  analysis  resulted  in 
backazimuth  estimates  correct  to  within  10  degrees  of  the  true  value. 
Note  that  in  most  cases,  the  very  early  portion  of  the  signal  jdelded 
the  best  estimate  of  the  backazimuth,  and  that  the  polarization 
characteristics  of  the  immediate  P-wave  coda  tend  to  fluctuate 
strongly. 
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Figure  16.  Time  series  data  for  Explosion  #9.  (Top)  Vertical  component, 
(middle)  radial,  (bottom)  transverse. 
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Figure  17.  Fourier  amplitude  spectrum.  Explosion  #9:  Vertical 
component  P  wave  arrival  (solid)  and  pre-P  wave  noise 
(dashed).  Ten  second  time  windows  were  used. 
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Figure  18.  Results  of  polarization  processing  for  Explosion  #9.  From  the 
top,  the  traces  are:  (1)  Vertical,  (2)  radial,  (3)  transverse,  (4) 
stadon-source  backazimuth  (deg),  (5)  recdlinearity,  and  (6) 
three>component  signal/noise  ratio.  Note  that  for  rotated 
components,  the  true  backazimuth  angle  is  180  degrees. 

Vertical  lines  indicate  time  at  which  the  recdlinearity  is 
maximum.  The  duration  of  the  moving  window  for  polarization 
analysis  is  0.5  seconds. 
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.  Time  series  data  for  Explosion  #12.  (Top)  Vertical 
(middle)  radial,  (bottom)  transverse. 
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Figure  20.  Fourier  amplitude  spectrum.  Explosion  #12;  Vertical 
component  P  wave  arrival  (solid)  and  pre-P  wave  noise 
(dashed).  Ten  second  time  windows  were  used. 
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Figure  22.  Time  series  dEta  for  Explosion  #17.  (Top)  Vertical  component, 
(middle)  radial,  (bottom)  transverse. 
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Figure  23.  Fourier  amplitude  spectrum.  Explosion  #17:  Vertical 
component  P  wave  arrival  (solid)  and  pre-P  wave  noise 
(dashed).  Ten  second  time  windows  were  used. 
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Figure  24.  Results  of  polarization  processing  for  Explosion  #17.  From  the 
top,  the  traces  are:  (1)  Vertical,  (2)  radial,  (3)  transverse,  (4) 
station-source  backazimuth  (deg),  (5)  rectilinearity,  and  (6) 
three-component  signal/noise  ratio.  Note  that  for  rotated 
components,  the  true  backazimuth  angle  is  180  degrees. 

Vertical  lines  indicate  time  at  which  the  rectilinearity  is 
maximum.  The  duration  of  the  moving  window  for  polarization 
analysis  is  0.5  seconds. 
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Figure  26.  Fourier  amplitude  spectrum.  Explosion  #33:  Vertical 
component  P  wave  arrival  (solid)  and  pre-P  wave  noise 
(dashed).  Ten  second  time  windows  were  used. 


32 


15-JUNE-92  23:26:20  UTC  RECT.  0.97:  BAZ.  184.7:  S/N  8.7 


Figure  27.  Results  of  polarization  processing  for  Explosion  #33.  From  the 
top,  the  traces  are:  (1)  Vertical,  (2)  radial,  (3)  transverse,  (4) 
station-source  backazimuth  (deg),  (5)  rectilinearity,  and  (6) 
three-component  signal/noise  ratio.  Note  that  for  rotated 
components,  the  true  backazimuth  angle  is  180  degrees. 

Vertical  lines  indicate  time  at  which  the  rectilinearity  is 
maximum.  The  duration  of  the  moving  window  for  polarization 
analysis  is  0.5  seconds. 


8.  Rg'Wave  Analysis 

The  complications  discussed  above  make  accurate  backazimuth 
determination  difficult  for  this  delay>fired  explosion  data  set.  In  an 
effort  to  extract  more  information  from  the  single-station  data,  we 
have  examined  the  potential  for  using  the  Rg  phase  to  determine 
source  direction. 

Several  of  the  explosion  signals  exhibit  Rg  phases.  At  BLA,  Rg 
is  a  Rayleigh  wave  with  group  velocities  in  the  range  1.8  to  2.3 
km/sec.  It  is  apparent  on  short  period  recordings  which  have  been 
bandpass  filtered  for  the  range  0.2  to  1.0  Hz. 

The  backazimuth  of  the  source  was  estimated  from  the  Rg- 
wave  by  noting  that  for  a  laterally  homogeneous,  isotropic  earth,  the 
phase  angles  of  the  Rg  radial  component  Fourier  harmonics  lead 
those  of  the  vertical  by  90  degrees.  In  the  presence  of  noise  and 
other  seismic  phases,  this  phase  difference  can  be  used  to  estimate 
the  backazimuth  if  energy  on  the  transverse  component  in  the  Rg 
frequency  band  is  sufficiently  incoherent  with  respect  to  that  on  the 
vertical  and  radial  components.  The  approach  we  have  used  is 
described  below. 

9.  Method  of  Rg-Wave  Polarization  Analysis 

Let  z(a)),  n(<a)  and  e(<a)  be  the  Fourier  transforms  of  the 
windowed  time  series  containing  the  Rg  arrival  on  the  vertical, 
north-south  and  east-west  components.  For  an  arbitrary  rotation 
angle  (>,  representing  a  trial  source-station  azimuth,  the  radial 
component  of  motion  is  given  by 

r(<o)  =  n(a))  cos(4>)  +  e(o))  sin  (<j>) ,  (3) 

and  the  phase  spectrum  is 

0^(0)  =  arg  r(a)).  (4) 
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Similarly,  for  the  vertical  component,  the  phase  is 

e  ^(o))  =  arg  z(o)),  (5) 

Because  the  Rg  signal  is  bandlimited,  we  define  oj  to  represent  the 
j'th  discrete  Fourier  frequency  in  the  range  con  ^  odj  ^  com.  j  =  n, 
n4-l....m,  where  con  and  (Om  are  the  lower  and  upper  bandlimits  of  the 
signal. 

Let  _ 

V2  2 

z(o>.)  +  r(ci)^)  ,  if  j  j^0j.(a)>0^(o)^)j-jc/2 1<  a 

=  0  ,  ifj  j^0j.(ti)>0^(a)^)j-ji/2|>a  (6) 

where  a  is  a  pre-selected  phase  difference  tolerance.  By  summing 
over  the  frequency  band  of  interest,  we  obtain  a  measure  of  the 
coherent  signal  amplitude  in  the  vertical-radial  plane.  This  sum  is 
normalized  by  the  total  three-component  amplitude  in  the  frequency 
band  to  yield 

m  y  [2  T~  2 

P=I  A.  /  /  z(a).)  +  n(to.)  +  e(o).)  .  (7) 

j=n  J  V  J  I  J  I  I  J 

This  quantity  is  unity  when  the  rotation  angle  ^  is  such  that  all 
horizontal  motion  is  in  the  radial  direction,  and  the  phase  differences 
between  radial  and  vertical  Fourier  harmonics  in  the  chosen 
frequency  band  are  in  the  range  ic/2  ±  a.  In  practice,  P  is  evaluated 
with  ^  ranging  from  0  to  359  degrees,  in  1  degree  increments,  with  a 
phase  difference  tolerance  of  10  degrees.  The  value  of  <)>  resulting  in 
a  maximum  value  of  P  is  the  estimated  Rg  azimuth  of  approach. 
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10.  Results  of  the  Rg-Wave  Analysis 

Twenty-seven  signals  exhibiting  Rg  arrivals  were  processed 
using  the  method  described  above.  The  calculated  P  values  ranged 
from  0.08  to  0.76.  Figures  28  through  31  show  the  results  of 
analysis  for  four  events.  Figures  28  and  29  show  results  for 
explosions  #22  and  #19,  respectively.  These  events  feature  the 
largest  amplitude  Rg  arrivals  in  the  data  set.  Figures  30  and  31  show 
results  for  explosions  #11  and  #29  which  feature  much  smaller  Rg 
arrivals.  Note  that  the  method  is  not  highly  dependent  upon  the 
amplitude  of  the  Rg  arrival.  This  is  further  illustrated  in  Figure  32 
which  plots  the  backazimuth  error  as  a  function  of  the  parameter  P. 
In  some  cases,  good  estimates  of  the  backazimuth  were  found  for  P 
values  as  low  as  0.08. 

The  backazimuths  estimated  from  the  Rg  phase  are  compared 
with  those  derived  from  the  P-wave  in  Figure  33.  Fot  Rg,  the  mean 
error  is  -2  degrees;  the  standard  deviation  is  26  degrees,  which  is 
only  S  degrees  larger  than  the  result  obtained  from  the  P-wave 
arrivals. 


8ACKA2IMUTH  (OEG) 


Figure  28.  Results  of  Rg  phase  analysis  for  Explosion  #22.  (Top)  Unfiltered  series 
(vertical,  radial,  and  transverse  components).  (Middle)  Bandpass  (0.2  to 
1.0  Hz)  filtered  time  series.  The  vertical  lines  indicate  time  segments  used 
for  analysis.  (Bottom)  Relative  in-phase  amplitude  (parameter  P  in 
Equation  7)  versus  backazimuth  (rotation  angle  #).  The  vertical  line 
indicates  the  true  backazimuth  of  180  degrees. 
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Figure  29.  Results  of  Rg  phase  analysis  for  Explosion  *19.  (Top)  Unfiltered  series 
(vertical,  radial,  and  transverse  components).  (Middle)  Bandpass  (0.2  to 
1.0  Hz)  filtered  time  series.  The  vertical  lines  indicate  time  segments  used 
for  analysis.  (Bottom)  Relative  in-phase  amplitude  (parameter  P  in 
Equation  7)  versus  backazimuth  (rotation  angle  *).  The  vertical  line 
indicates  the  true  backazimuth  of  180  degrees. 
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Figure  31.  Results  of  Rg  phase  analjrsis  for  Explosion  #29.  (Top)  Unfiltered  series 

(vertical,  radial,  and  transverse  components).  (Mid^e)  Bandpass  (0.2  to 
1.0  Hz)  filteied  time  series.  The  vertical  lines  indicate  time  segments  used 
for  analysis.  (Bottom)  Relative  in-phase  amplitude  (parameter  P  in 
Equation  7)  versus  backazimuth  (roution  angle  9).  The  vertical  line 
indicates  tte  true  backazimuth  of  180  degrees. 
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Figure  32.  Rg  phase  backazimuth  error  versus  relative  amplitude  parameter  P. 
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Figure  33.  Backazimuth  errors  derived  from  P  and  Rg  phases. 
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11.  Conclusions  and  Recommendations 

We  have  compared  the  single-station  backazimuth  estimates  at 
BLA  with  independently  derived  locations  for  a  set  of  mining 
explosions  and  find  that  for  sources  in  the  northwest  quadrant, 
polarization  analysis  of  the  P-wave  arrival  gives  an  essentially 
unbiased  estimate  of  the  source  direction  (mean  error  6  degrees),  but 
that  the  scatter  in  the  results  is  substantial  with  a  standard  deviation 
for  35  measurements  of  21  degrees. 

Characteristics  of  the  P-wave  signal  which  contribute  to  the 
errors  in  the  estimates  include:  (1)  emergent  initial  motions  from  the 
delay-fired  explosions;  (2)  off-azimuth  arrivals  and/or  converted 
phases  arriving  early  in  the  P  coda;  (3)  steep  angles  of  incidence, 
averaging  22  degrees  from  vertical,  which  tend  to  reduce  the 
signal/noise  ratios  on  the  horizontal  components.  Optimum  data 
segments  for  reliable  backazimuth  estimates  are  restricted  to  short 
time  intervals  (1  second  or  less)  beginning  with  the  initial  P-wave 
motion. 

The  Rg  phase  is  present  in  many  of  the  explosion  signals  at 
BLA.  Analysis  of  this  arrival  using  a  phase  difference  criterion 
results  in  backazimuth  estimates  which,  on  average,  have 
approximately  the  same  reliability  as  those  derived  from  the  initial 
P-wave  motion.  It  appears  that  this  approach  can  often  provide  an 
alternative  means  of  estimating  the  source  direction,  which  could  be 
important  in  situations  where  the  P-wave  signal/noise  ratio  is  small. 
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